Supersolids in one dimensional Bose Fermi mixtures 
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Using quantum Monte Carlo simulations, we study a mixture of bosons and fermions loaded on 
an optical lattice. With simple on-site repulsive interactions, this system can be driven into a solid 
phase. We dope this phase and, in analogy with pure bosonic systems, identify the conditions under 
which the bosons enter a supersolid phase, i.e., exhibiting at the same time charge density wave and 
superfluid order. We perform finite size scaling analysis to confirm the presence of a supersolid phase 
and discuss its properties, showing that it is a collective phase that also involve phase coherence of 
the fermions. 
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The supersolid (SS) phase in which superfluid and solid 
order coexist was first examined in He systems more 
than 50 years ago^. The presence of a SS phase was con- 
clusively demonstrated in several lattice models such as 
variants of the bosonic Hubbard model. However, this 
phase has not been unambiguously observed experimen- 
tally and, with the recent progress in low temperature 
physics, the search for experimental evidence of a SS 
phase has been reinvigorated. 

In solid helium, a non-classical reduction of the mo- 
ment of inertia was observed in torsional oscillator 
measurements This reduction is due to the appear- 
ance of a superfluid fraction in the material and has 
been first interpreted as a sign of a SS phase. However, 
several experimental^!^ and numerical studiesi^^ sug- 
gested that this superfluid behavior was in fact due to the 
presence of non local defects in the system (grain bound- 
aries for examples) along which superfluid currents exist. 
The presence of a bulk supersolid phase thus appears to 
be questionable in these systems. 

Another promising approach for finding a SS phase 
emerged in the context of cold atoms loaded on optical 
lattices. Such systems, being well described by bosonic 
Hubbard models, are a good starting point to look for 
supersolids. However the conditions necessary for super- 
solids to appear in these models are not easily achieved in 
real experiments: One generally observes a direct solid- 
superfluid transition^ or a coexistence of superfluid and 
solid (in the case of a first order phase transition)^. Su- 
persolids are stabilized in these models by specific long 
range interaction a 12 i 13 ' 14 i 15 ' 16 i 17 i 18 or long range hopping 
termsJ^ Engineering precisely the values of the interac- 
tion as a function of distance appears to be tricky to 
achieve. 

Recently Bosc-Fcrmi mixtures have been introduced 
as a way to study the physics of fermions, us- 
ing sympathetic cooling with the bosons to reach 
very low temperatures! 20 ! 21 However several theoreti- 
cal studie9&22i2£i2&%L21 have shown that such mixtures 
loaded on optical lattices have a rich phase diagram 
where collective phases of fermions and bosons appear. 



Interestingly, with simple repulsive on-site interactions 
between fermions and bosons, the system can be driven 
in a solid phase where density- wave order develops^ 
We will study here the doping of such a solid phase 
and see under what conditions it could be driven into 
a bosonic supersolid phase. Other recent studies sug- 
gested the presence of a SS behavior in Bose-Fermi 
mixtures 2§£&2ML 

The paper is organised as follow. In section I, we in- 
troduce the model and discuss the solid phase which is 
the starting point of our study. In section II, we explore 
the different ways of doping this system and determine 
which one could lead to a supersolid phase. Finally, in 
section III, we perform finite size scaling analysis of phys- 
ical quantities to verify if the supersolid phase persists in 
the thermodynamic limit. 



I. SOLID PHASE IN THE BOSE-FERMI 
HUBBARD MODEL 

We study a one dimensional Hubbard model for a mix- 
ture of bosons and polarized (spinless) fermions. The 
Hamiltonian is given by: 



H = f~*b b t+l b r - *f /r+l/r + b.C 



+ £ \U hb nhr{n \ l) + U b{ ny^J (I) 

where b\ (b r ) creates (destroys) a boson on site r while 
/t and f r are the corresponding operators for fermions. 
The first term in Eq.JT]) describes tunneling of bosons 
and fermions between neighboring sites, with different 
associated energies th and U . In the following, the energy 
scale is fixed by choosing t\> = I. nj? and n r are the 
bosonic and fermionic number operators at site r and Uhb 
and C/bf are the boson-boson and boson-fermion contact 
repulsion terms. L is the total number of sites in a one- 
dimensional chain. 
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We will be mostly interested in the behavior of bosons 
and will, therefore, study the bosonic Green function 
Gb(-R) which measures phase correlations: 

1 L 

G h (R) = I J2( b t + R b r + bib r+R ) (2) 
r=l 

as well as the boson superfluid density p s = 
(W 2 L)/(2f3tb) where W is the winding number— and /3 
the inverse temperature. The density-density correlation 
function 

MR) = K n r + R n r) ( n r+n) <"r>] (3) 

r=l 

and its Fourier transform, the structure factor Sh(k), give 
information on possible solid density-wave order. 

We study this model using two different versions of 
the worm algorithm: the canonical worm algorithm 
(CW) 33 ' 34 and the directed stochastic Green function al- 
gorithm (DSGF)^^. The CW method, which is very 
efficient for measurement of equal time Green functions, 
has to be modified to include simultaneous updates of 
bosons and fermions to allow the study of mixtures 3 - 7 -. 
However, this algorithm sometimes becomes inefficient 
for large system sizes (L > 20) especially failing to sam- 
ple different winding numbers which is necessary to cal- 
culate the superfluid density 3 -^. In those cases we used 
the DSFG which explores a much larger configuration 
space and thus allows more efficient fluctuations of the 
winding number and the measure of p s . The two algo- 
rithms gave results in agreement for other quantities in 
the range of sizes we used. 

This model was widely studied in the special case 
where the total number of particles equals the number 
of sites JVb + Nf = L where iV"b(f) is the number of 
bosons (fermions) i 24 i 37 For large enough values of the re- 
pulsions, double occupancy of sites is forbidden and one 
can describe the system in terms of a pseudo spin-1/2 
a* = — n\. = ±L 22 i 23 The Hamiltonian can then be 
mapped, in the low energy limit, into an effective Heisen- 
bcrg Hamiltonian 

H eS = E J *y + + E J *«+i (4) 

r r 

where J xy = -t h t { /U bi and J z = (ig + tj)/{2U hi ) - 
t^/Uhb- When J z > \J xy \, the pseudo-spin system en- 
ters a Neel antiferromagnetic phase along the z axis^ 
In terms of bosons, this antiferromagnetic order corre- 
sponds to a density wave order with alternating occupied 
and empty sites (see Fig. [I]). The Green function, G^{R), 
decays exponentially indicating the absence of phase co- 
herence. This phase persists in the thermodynamic limit 
at zero temperature (see Fig. [1] inset). 
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FIG. 1: (Color online) Density-density correlation D^{R) 
(multiplied by ( — l) fl )and Green function Gb(R) for the 
bosons in the solid ordered phase. Inset: Finite size behavior 
of the bosonic structure factor Sb(w) for the same parameters 
and sizes ranging from L — 10 to 70. 



II. DOPING THE SOLID PHASE 

A similar density wave ordered phase for N b = L/2 is 
observed in a one-dimensional bosonic system with near- 
neighbour repulsion^ In this latter system a supersolid 
phase can be present when the solid phase is doped by 
adding bosons. However, for a supersolid phase to ap- 
pear, the interactions must be chosen so that the added 
bosons do not introduce defects in the previous solid or- 
der: They must preferentially come on sites already oc- 
cupied by bosons. A mean field argument yields that the 
repulsion between bosons located on neighboring sites 
must be greater than half the on-site repulsion Uhb- 

Following this example, we consider a system where 
the repulsion between bosons, t/bb = 15, is smaller than 
the boson-fermion repulsion Uhf = 18. Starting from the 
solid phase obtained for Ab = Nf = L/2 and if = 4^ 
we changed slightly the number of the different types of 
particles and observed how the solid order was modified 
(see Fig. 0. 

We found that upon changing the number of fermions, 
some defects are introduced in the solid order, which 
are exposed by the characteristic beating in the density- 
density correlations (Fig. [21 top)^ On the other hand, 
when one changes the number of bosons, the density wave 
order of alternating empty and filled sites persists (Fig. [21 
bottom). A surprising result is that this wave order is 
present even when the number of bosons is reduced be- 
low half filling, unlike what happens in the purely bosonic 
model with near neighbor repulsioiii 17 

Examining the bosonic structure factor 5b (n) (Fig. [3]), 
we observe that only the case Nf = L/2 leads to large 
Shin) and therefore to the long range density order neces- 
sary for the establishment of SS. We also observe in Fig. [3] 
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FIG. 2: (Color online) Density-density correlations for differ- 
ent boson fillings, iVb, and for two different fermion fillings, 
Nf — L/2 — 2 and N{ = L/2. In the first case (a), we observe 
the beating characteristic of the pseudo spin phase at non zero 
magnetisation along z with a maximum when Nf + Nb = L. 
In the second case (b), we see a long range order (or quasi 
long range) for the bosons. 



that when the bosonic population is doped above or be- 
low half filling, Sb(ir) drops but remains rather apprecia- 
ble especially above half filling. Finite size scaling is re- 
quired to establish if these non- vanishing values of Sb(7r) 
persist in the thermodynamic limit (see section UIT|) . 

As expected, quasi long range phase coherence is recov- 
ered as soon as the solid is doped away from half-filling 
(see FigU]) and the bosons become superfiuid. This phase 
coherence is stronger when the system is doped above 
half-filling. 

These results lead to the conclusion that a supersolid 
phase can only be found for Nf = L/2, where no defects 
are introduced in the solid order. Doping the bosons 
above half-filling provides the best candidate to observe 
a supersolid as both the structure factor and the phase 



FIG. 3: (Color online) Structure factor Sb(?r) as a function 
of the number of bosons for different fermion fillings, Nf. 
The oscillations characterizing a density order only develop 
for N f = L/2. 




FIG. 4: (Color online) Phase correlation, Gh(R), for differ- 
ent boson fillings. There is always an algebraic decay char- 
acteristic of a superfiuid in one dimension, except at double 
half-filling in the solid phase. 



coherence are larger in this case but doping below half- 
filling could also yield to a supersolid phase. Finite size 
analysis of the behavior of p s and Sb(ir) is needed to ver- 
ify if the supersolid phase persists in the thermodynamic 
limit. 



III. FINITE SIZE ANALYSIS 

In this system, the sizes that can be used in the finite 
size scaling analysis arc quite limited, fn order to avoid 
sign problems, Nf = L/2 must be odd. In addition, L 
cannot be very large in order for the simulation to con- 
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FIG. 5: (Color online) Finite size scaling of the structure 
factor Sb(7r) for different values of the hopping parameter if. 
The structure factor goes to zero in the large size limit for 
tf = 4. For if > 4, in the range of sizes that are accessible, 
the structure factor is non zero. 



verge in reasonable time. We used two different densities 
of bosons, Nb/L = 0.4 and 0.6 (one below and one above 
half-filling) which give an integer number of bosons for 
sizes L — 10, 30, 50 ■ ■ ■ When we could not obtain ex- 
actly these densities with the given constraints, we used 
the number of bosons that gives the closest density. 

We first performed the finite size analysis with the pa- 
rameters used in the first part of this paper (Uhb = 18, 
U hf = 15, if = 4) and above half-filling (N h /L = 0.6). 
Figures [5] and [6] show that, for this case, the structure fac- 
tor goes to zero in the thermodynamic limit while the su- 
pcrfluid density remains finite. This indicates that what 
appears to be a supersolid phase for small L is in fact 
a superfluid. Varying the different available parameters 
(Ubb, Uhf, if), we observed that increasing if increases 
noticeably the value of Shift), whereas varying the in- 
teraction did not yield similar variations. For if > 5 
(see Fig. [5]) Sb(ft) extrapolates to a finite value as L 
increases. This indicates that for these values of if the 
density wave order survives in the thermodynamic limit 
when iV b = 0.6L and 7V £ = L/2. 

In Fig. [6] we see that the superfluid density changes 
very little with L or if: p s then goes to a finite value 
when L — > oo. As mcntionned above, obtaining precise 
values for p s is difficult but, on the other hand, the Green 
function, Gb{R), is measured with very good statistical 
accuracy and can, therefore, also be used to character- 
ize the nature of the phase coherence. Figure [7] shows 
that Gb(R) exhibits power law decay for all the sizes and 
values of the fermion hopping parameter we have stud- 
ied (for A^b = 0.6L), confirming the existence of a quasi 
condensate which leads to a superfluid behavior in the 
presence of long range density order for the fermions. 

This indicates the presence of a supersolid phase in this 
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FIG. 6: (Color online) Finite size scaling of the superfluid 
density p a in the case where the system is doped above half- 
filling. Reliable results for this quantity are difficult to obtain 
for sizes larger than L > 42. However, p s shows very little 
variations when L or if vary and should then remain non-zero 
in the large L limit for all cases. 



So far, we have concentrated on the behavior of the 
bosons. However, the supersolid phase we have exposed 
is necessarily a collective phase of fermions and bosons. 
To study the behavior of fermions or the collective be- 
havior of fermions and bosons, we introduce the Green 
functions for fermions Gf, boson-fermion pairs Gbf, and 
boson- hole pairs Gbh- 



L=30, 50, 70 (or 62), 



=L/2, N b =0.6L 




FIG. 7: (Color online) Finite size scaling of the Green function 
Gb(-R) for several values of the hopping parameter if. Gb(-R) 
always appears to decay algebraically. The curves have been 
multiplied by a factor indicated in the legend between paren- 
theses. 
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FIG. 8: (Color online) Green functions for bosons, fermions, 
boson-fermion pairs, and boson-hole pairs in the candidate 
supersolid phase. All the Green functions have algebraic de- 
cay. 
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FIG. 9: (Color online) Finite size scaling of the structure 
factor SbM below half-filling (N b /L = 0.4) for tt = 7. The 
structure factor goes to a non zero value in the large size 
limit. The observed oscillations are due to the fact that we 
obtained exactly the desired density (Nb/L = 0.4) only for 
L = 30, 50, 70. 



1 

Gi(R) = ^(fr+Rfr + fifr+R 
r=l 

1 L 

G h{ (R) = ^(/^U/' + L 

r=l 
1 L 



(5) 
(6) 
(7) 



Figure [5] shows that these Green functions have algc 



braic decay in the candidate supersolid phase. The dom- 
inant correlations are those of individual bosons, they 
have the slowest decay and the largest values. However, a 
description only in terms of bosons does not fully charac- 
terize this phase; we see that individual fcrmion, boson- 
hole and boson-fermion pairs also have algebraic phase 
correlations and are also relevant degrees of freedom in 
the system. The boson-hole pairs have much stronger 
phase correlations than the boson-fermion pairs. This is 
to be expected since anticorrelated movement of bosons 
and fermions is a generic phenomenon in such mixtures 
with repulsive interaction a 22 ' 25 ' 40 . 

In a similar way to the case above half-filling, we per- 
formed finite size scaling analysis when the system is 
doped below half- filling (Nb/L = 0.4). We studied a case 
where the value of if is quite large (U = 7) since we ob- 
served in the previous case that it improves the structure 
factor. We find that, as for the case above half filling, 
the structure factor and the superfluid density both go 
to a finite value in the large L limit and that a supersolid 
phase is, therefore, thermodynamically stable (see Fig. 
[H]). This behavior is different from what happens in the 
bosonic case with near neighbor interactions where, dop- 
ing below half filling leads to a superfluid with soliton-likc 
quasiparticlesi 17 i 39 



IV. CONCLUSION 

In this paper we have studied, using exact QMC simu- 
lations, the formation of supersolid phases in the ground 
state of Bose-Fermi mixtures on optical lattices. At dou- 
ble half filling, iVb = N[ = L/2, and for sufficiently 
large tf/th, the system exhibits long range density or- 
der as exposed by Sh(Tr)~ Inspired by the behavior of 
the extended bosonic Hubbard model, which exhibits a 
supersolid phase above half filling when the near neigh- 
bor repulsion is large enough compared to the contact 
term ; 12 ' 17 , we found that, similarly, when Ubi > U b b and 
the system is doped by adding bosons, the system en- 
ters a supersolid phase. It is important to keep in mind 
that this phase is a collective Bose-Fermi phase: the 
various Green functions (boson-boson, fermion-fermion 
and boson-fermion) we have presented demonstrate this 
clearly. 

Surprisingly, we have also found that upon doping the 
system below half filling for the bosons, we also obtain 
a stable supersolid phase. This is different from the be- 
havior of the purely bosonic system. 
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